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Abstract 

C\ , • We consider the one-dimensional porous medium equation u t — {u n u x ) x + ^-u n u x . We 

derive point transformations of a general class that map this equation into itself or into 
equations of a similar class. In some cases this porous medium equation is connected 

rp with well known equations. With the introduction of a new dependent variable this 

partial differential equation can be equivalently written as a system of two equations. 
Point transformations are also sought for this auxiliary system. It turns out that 
in addition to the continuous point transformations that may be derived by Lie's 
method, a number of discrete transformations are obtained. In some cases the point 
transformations which are presented here for the single equation and for the auxiliary 

^^ \ system form cyclic groups of finite order. 
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1 Introduction 



Probably the most useful point transformations of partial differential equations (pdes) 
are those which form a continuous (Lie) group of transformations, each member of which 
leaves an equation invariant. In addition to possessing continuous groups of symmetries, 
many pdes also possess discrete symmetries. However, the classical method of deriving the 
continuous symmetries, which is based on the linear nature of infinitesimal transforma- 
tions, is not directly applicable to the derivation of groups of discrete symmetries. Such 
S^ example is given by Kingston and Sophocleous [1] who found that the reciprocal point 

transformation (double application gives the identity transformation) x' = x/t, t' = 1/t, 
v! = —ut + x leaves the Burger-type equation u t + uu x + (f(t) — f(l/t))u xx = invariant, 
which is a symmetry additional to the Lie point symmetries obtained from the classical 
approach. 

Point transformations may also be used to link a pde with an equation of a different 
form. For example, it is sometimes possible to connect an equation to a canonical form 
for which there is an established theory. Such cases are given in the present paper. 

For these reasons it is advantageous to study point transformations directly in finite 
forms with the ultimate dual aims of finding the complete set of point transformation 
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symmetries of pdes and also discovering new links between different equations. A good 
example, of the approach which is used here, is the work by Sophocleous and Kingston who 
classify the point transformations of four common classes of one-dimensional non-linear 
wave equations [2]. 

We consider the one-dimensional porous medium equation of the form 

u t = (u n u x ) x + -u n u x . (1.1) 

ix x 

We introduce point transformations of a general class and investigate when such transfor- 
mations map Eq.(l.l) into itself. In addition we classify point transformations that link 
Eq.(l.l) to different equations but of the same class, where in some cases it turns out to 
be well known equations. A number of the point transformations obtained form cyclic 
groups of finite order. 

The next section is devoted to the above analysis, where firstly it makes use of a result 
obtained by Kingston and Sophocleous [3] which is summarised in the following theorem: 

Theorem. The point transformation x' = P(x,t,u), t! = Q(x,t,u), v! = R(x,t,u) 
transforms the pde 

u' t , = H'(x',t',u',Ui, . . . ,u' n ) 

to the pde 

u t = H(x,t,u,u 1 ,...,u n ) 

where u n = -q^c, u' n = ^zik, n > 2 and H and H' are polynomials (non-negative integral 
powers) in u±, . . . ,u n and u' 1: . . . ,u' n respectively, if and only if 

P = P(x,t), Q = Q(t) 

and the functions H and H' satisfy the relation 

PxRuH = P x QtH + PtR x + PtRuUx — PxRt- (1-2) 

The proof of this theorem can be found in ref. [3] . 

Equation (1.1) can be written in a conserved form and therefore with an introduction 
of a new dependent variable may be represented by a system of two pdes. In the Sections 3 
and 4 we search for point transformations of a general class that map this system into 
itself or to a system of the same class. In some cases, such point transformations of this 
auxiliary system lead to non-local transformations of Eq.(l.l). 

Equation (1.1) is a special case of the generalised porous medium equation 

u t = (u n u x ) x + f(x)u s u x + g{x)u m . (1.3) 

The classical point symmetries have been presented by Gandarias [4]. Further study, along 
the lines of this paper, of Eq.(1.3) may therefore be useful. 
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2 Point transformations 

We consider point transformations of the form 

x' = P(x, t, u), t' = Q(x, t, u), u' = R(x, t, u) (2.1) 

which transform the pde 

u t> = (/(w')«x')x' + 9 ( x 'i u '^' ( 2 - 2 ) 

to the pde 

Ut = (F(u)u x ) x + G(x,u)u x , (2.3) 

where F not a constant. We assume that point transformation (2.1) is non-degenerate in 
the sense that 

d(P, Q,R) ^ n d(P(x,t,u(x,t),Q{x,t,u{x,t))) ^ n 

-— — ^ and — 7^ 0. 

d{x,t,u) r d(x,t) r 

Firstly, we note that each linear transformation of the form 

x' = cix + c 2 , t' = c 3 t + c 4 , v! = c 5 u + c 6 , (2.4) 

where q are constants (possibly complex), C1C3C5 / 0, achieves this purpose. That is, 
Eq.(2.2) is transformed to Eq.(2.3) with the function F(u) defined in terms of f(u') and 
G(x, u) in terms of g(x', v!) by 

F(u) = c 3 q 2 /(c 5 n + c 6 ), G(x,u) =c 3 c^ 1 g(cix + C2,c 5 u + c 6 ), (2.5) 

respectively. Hence, point transformations of the form (2.4) may be filtered out of equa- 
tions arising in the discussion, without loss of generality, with the understanding that all 
point transformations obtained may be augmented by Eqs.(2.4) and (2.5). 

We now turn into the special case of Eq.(2.3), where F(u) = u n and G(x,u) = -u n , 

u t = {u n u x ) x + ^u n u x . (2.6) 

The Lie point symmetries of the equation u t = (u n ) xx + ^ {u n ) x , which is equivalent to 
Eq.(2.6) with the exception of the case where n = — 1, have been classified by Gandarias 
et al [5]. The symmetries of Eq.(2.6) are summarised in the following table, where (X, T, U) 
are the components of the symmetry generator X4^ +T-q^ + [/^: 





n 
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We search for point transformations of the form (2.1) that map Eq.(2.2) into Eq.(2.6). 
We note that both of these equations are such that the Theorem, given in the Introduction, 
can be applied. Therefore we have Q = Q(t) and P = P(x,t). Also we substitute 
H = u n u xx + nu n - l u 2 x + %u n u x , H' = f(u')u x , x , + ^u' x , 2 + g(x', u')u' x , into Eq.(1.2) and 
with the use of the expressions of u' x ,, u' x , x , [3] the relation (1.2) becomes 

xuP x R u (u n P 2 - Q t f) u xx 

+xP x (nu n P 2 R u - uQ t R uu f - u Q t R 2 u ^-\ u\ + u[xQ t P xx R u f - xP t P 2 R u 



dA 

du'J 
+xu ( Q t P xx R x f - PtP A x R x + P 6 x Rt - QtP x R x g - QtPxRxxf - QtPxR 



+fiu n P^R u - xQ t P 2 R u g - 2xQ t P x R xu f - 2xQ t P x R x R u ^)u 



(2i 



du' 

We set, successively, the coefficients of u xx , u x , u x and the term independent of u x and 
u xx in Eq.(2.8) equal to zero to obtain four identities. These identities enable the desired 
point transformations to be derived and ultimately impose restrictions on the functional 
forms of Q(t), P(x,t), R(x,t,u), f(u') and g(x',u'). For example, from the first two we 
deduce that / = -^u n and also that R is linear in u. From these we conclude that /(it') 
is linear in u' n , but using the transformation (2.4) we can simply write /(it') = it" 1 . 

In the following analysis we omit any further calculations. It turns out that this analysis 
can be split into two cases:(l) n = — 1; (2) n ^ — 1. Firstly, we point out that the discrete 
transformations which correspond to the symmetries I 1 ]— T3 admitted by Eq.(2.6) are of 
the form (2.4). 

Case 1: (n = — 1) Equation (2.6) takes the form 

Ut =(^) +JL Ux . (2.9) 

V U / x xu 
If n = 1, Eq.(2.9) is mapped into itself by the transformation 

x' = x k+1 , t' = t, u' = - l , x- 2k u. (2.10) 

This latter transformation is the continuous symmetry which is represented by the sym- 
metry T5 = xlnx^ + 2(lnx — l)itg- admitted by Eq.(2.9). 
If we apply transformation (2.10) A^-times we obtain 

X (N) = x(fc+ l)* t (N) = t U (N) = 1 -2( fc+ l)"+2 U- (2n) 

[k + l) /iV 

We note that this point transformation forms a cyclic group of order N if k is a root of 
the equation 

(k + 1) N - 1 = 0. (2.12) 

For example, if N = 2 then k = —2 (we omit the obvious solution k = 0) and we get the 
transformation x' = 1/x, t' = t, u' = x 4 u which forms a cyclic group of order 2. In fact, it 
is trivial to note that k = and k = —2 are the only real roots of equation (2.12). If we 
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introduce complex roots of Eq.(2.12), then we may construct cyclic groups of any finite 
order. For example, the point transformation x' = x 1 , t' = t, v! = —x 2 ~ 2l u forms a cyclic 
group of order 4. 

Furthermore identity (2.8) produces the following three equations that mapped into 
Eq.(2.9) (fj, = 1) by the corresponding discrete transformations: 

v '=(!!^.) ; x' = Hn|x|, t' = t, v! = — U X 2 \ (2.13) 

\ u 1 „, ' k z 



11, 



i \ / 1 i \ „k 



'-I^A xl-Lxl-r'l,,'. -r'-fL *'-W „' - _L„r2(l-k) . 



■f)A^ + 2 x ') u '"'' x ' = Vf l ' = lnt - u ' = v u * i(2 ' 14) 

/ _ / U a;' 1 ,1 /„,/ . , _ fc In |x| , ,1 2 



?i t / 



, H — x'u'r, x' = 7= — , t' = hit, u = -r^ux . (2-15) 

w J x , 2 y/t k z 



We note from result (2.13) that Eq.(2.9) is connected with a well known nonlinear dif- 
fusion equation. Therefore any solution of this diffusion equation can be transformed into a 
solution of Eq.(2.9) using (2.13). For example, point transformation (2.13) maps the solu- 
tion v! = 2t' (x' 2 + At' 2 \ of u' t , = (u' x ,/u') x > into the solution u = 2t [x 2 (ln 2 \x\ + At 2 )] ^ 
ofEq. (2.9). 

Finally, if \i = 0, the following equation with the corresponding point transformation 

<> = {^) x +(-^, + r)<>i *' = ^ f = ln*, u' = ue- 2 * (2.16) 

is mapped into Eq.(2.9) which has the form of the nonlinear diffusion equation ut = 
(u x /u) x . Result (2.16) can also be obtained by combining the results (2.13) and (2.14). 
Case 2: (n ^ —1) The continuous symmetry which corresponds to the infinitesimal 

n o 2(n+l) 

symmetry F^ = (n + 2)x n + 2 5 lux < n + 2 ) 5- admitted by Eq.(2.6) is the following 

x' = I x ™+2 + Cj , t = t, u ,=ux™+ 2 (x™+ 2 +C) , (2.17) 

where also n / —2, and \x = 3 ^+ 2 4 . 

In addition to the above result, identity (2.8) produces three more discrete transforma- 
tions. When \x ^ — ^— , the point transformation 

2n + 2 M^+n+2 M-i 

x = x 2™+ 2 , t =t, u =x"+ 1 u (2-18) 

/in + n + 2 

maps the equation 

m 

it 

<' = { u ' n <>) x/ + X — <> (2-19) 

into Eq.(2.6), where A = ~"^ 2 ^ . If A = \i = — ^ then the discrete point transfor- 

1 4 

mation x = -,£' = t, u' = x~™u, which forms a cyclic group of order 2, maps Eq.(2.6) 
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into itself. We observe that ft = ^+2 implies A = and Eq.(2.19) becomes the stan- 
dard nonlinear diffusion equation which is transformed to Eq.(2.6) (n 7^ — 2, —1) by the 
transformation 

71 + 2 2n + 2 _^_ 

x = x n + 2 , t =t, u = x n + 2 u. 

2n + 2 

Also fx = implies A = ^" 2 and we obtain the inverse of this latter point transformation 
which maps Eq.(2.19) into the nonlinear diffusion equation. 
When \x = — ^— , the point transformation 

x' = ln\x\, t' = t, u ' = x~n U (2.20) 

maps the equation 

u t' = \ u u x') x > + 2 —^- u u x> I 2 - 21 ) 

into Eq.(2.6). Eq.(2.21) is the Boussinesq equation of hydrology involved in various fields 
of petroleum technology and ground water hydrology. We observe that when n = —2 
(/i = 0) Eq.(2.21) is mapped by (2.20) into the well known nonlinear diffusion equation 
ut = {u~ 2 u x } . This diffusion equation has a number of properties. For example, admits 
Lie-Backlund transformations and there exists a transformation that maps it into the 
linear heat equation ut = u xx [6]. Also we point out that if n = — 1 we obtain the result 
(2.13) of the case 1. 

Finally, when \i 7^ — ^-, the point transformation 

fin-\-n-\-2 

, 2n + 2 x 2 ™+ 2 . M-i . 

x' = ^, t' = lnt, u' = x"+iu (2.22) 

fin + n + 2 -y/t 

maps the equation 

A^r + ^A <' (2-23) 

into Eq.(2.6), where A = n nn+n+2 ^ ' ■ -^ n ^ ne case wnere A* = ~~ ^~i identity (2.8) gives 
n = — 1 (fj, = 1, A = 0) and we simply reproduce the result (2.15). 

3 Potential transformations 

If we introduce the potential v, Eq.(2.6) can be written as a system of two pdes, 

v x = ux, v t = xu n u x -\ — -u n+1 , (3.1) 

n + 1 

if n 7^ — 1 and as 

x 
v x = ux, vt = —u x + (u — 1) In \u\, (3-2) 

u 

if n = — 1. 
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Similarly, if we write 

gixW)= m±^i, ( 3.3) 

X 

where h'(u') = -7-7, then Eq.(2.2) can also be written as a system of two pdes 

v' x , = x'u', v' t , = x'f(u')u' x , + h(u'). (3.4) 

We also consider the point transformation 

x' = P(x, t, u, v), t' = Q(x, t, u, v), u' = R(x, t, u, v), v' = S(x, t, u, v) (3.5) 

relating x, t, u(x,t), v(x,t) and x', t', u'(x',t'), v'{x',t') and we assume that it is non- 
degenerate. 

It is well known that in many cases pdes which can be written in a conserved form 
admit nonlocal symmetries, known as potential symmetries [7]. For example, in order to 
find potential symmetries for Eq.(2.6) we search for point symmetries for the system (3.1) 
(and (3.2)). If at least one of the infinitesimal generators of the variable x, t, u depends 
on the potential variable v, then the point symmetry of the auxiliary system (3.1) is also 
a potential symmetry of the Eq.(2.6). Otherwise, the point symmetry of (3.1) projects 
onto a point symmetry of Eq.(2.6). 

The symmetry analysis of the system (3.1) was carried out in the ref. [5]. It was shown 
that only one symmetry of (3.1) induces a potential symmetry admitted by Eq.(2.6). The 
rest of the symmetries project onto point symmetries of Eq.(2.6). The potential symmetry 
occurs if n = —2 and \x = — ^ and is given by 

2xv- 2{x 2 u 2 + uv)— + v 2 — . (3.6) 

OX OIL OV 

Symmetry (3.6) generates the one-parameter continuous group of transformations 

, x , , u(l — ev) 3 , v . _. 

x' = - ^, t' = t, u'= 2 V ' , v' = - . 3.7 

(1 — evy 2ex z u + 1 — ev 1 — ev 

Our goal in this section is to derive point transformations of the class (3.5) which map 
Eqs.(3.4) into Eqs.(3.1) (and Eqs.(3.4) into Eqs.(3.2)). Using (3.5) and the forms of u' x ,, 
u' t ,, v' x ,, v' t , in terms of u x , ut, v x , vt and derivatives of P, Q, R, S [8] and eliminating v x 
and vt from Eqs.(3.1) (or Eqs.(3.2)), Eqs.(3.4) become two identities of the form 

Ei(x,t,u,v,u x ,u t ) =0, E 2 (x,t,u,v,u x ,u t ) = 0, (3.8) 

where x, t, u, v, u x and ut are regarded as independent variables and E±, £2 are explicit 
polynomials in u x and Ut- 

Now, the coefficients of u x , u x , u% in E\ = and u x , ut in E 2 = 0, which must be 
identically equal to zero, give P u = Q x = Q u = Q v = S u = 0. Therefore, the point 
transformation (3.5) can be written in the simplified form 

x' = P(x,t,v), t' = Q(t), v! = R(x,t,u,v), v' = S(x,t,v). (3-9) 

To ensure that the point transformations are non-degenerate we require to have 

Q t R u (P x S v - P V S X ) + 0. (3.10) 
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In addition, we need to have 

P 2 + R 2 v ^0 (3.11) 

because otherwise the point transformations derived are equivalent to point transforma- 
tions which connect Eq.(2.2) and Eq.(2.6). Transformations that satisfy condition (3.11) 
shall be called potential transformations. 

The restricted forms (3.9) simplify the identities (3.8). For the convenience of the 
reader we omit any further calculations. We present the three potential transformations 
obtained: (1) n = -1; (2) n = -2; (3) n / -2, -1, 0. 

Case 1: If n = — 1 then the point transformation 



t v 4 ee -( ,; + 2 *) 
x' = e t+ 2, t' = t, u = - , v' = 2(elnx-t), (3.12) 



transforms the pdes 

x 1 
v' x , = x'u', v' t , = e—u' x , + 2(e - 1), (3.13) 

where e is a constant, to the pdes (3.2) where fj, = 1. Application of the point transfor- 
mation (3.12) 2N-times gives 

x , { 2N) = ^ t>{2N) = ^ u , {2N) = ffljMi v , { 2N) = e M y + ^M _ ^ (314) 

Therefore transformation (3.12) forms a cyclic group of order 2N if e is a root of the 
equation e = 1. Clearly, the only two real roots are e = 1 and e = — 1 where in these 
cases transformation (3.12) forms a cyclic group of order 2 and order 4, respectively. If 
we allow e to be complex, then we may construct cyclic groups of any even order. For 
example, if e = i then the point transformation x' = e t+v ' 2 , t' = t, u' 



ii 



ux 2 cxp(v+2t) ' 

v' = 2(ilnx — t) forms a cyclic group of order 8. 

If e = 1 the reciprocal transformation (3.12) maps Eqs.(3.2) (/j, = 1) into itself and 
therefore we also observe that the point transformation 



.r 



' = J+v/2 ,1 = . i 



t' = t, v' = 2(lnx-t) 



leaves the integrated form of Eq.(2.9) {u 



x 
vt = —v xx - 1 

Vx 

invariant. 

Case 2: If n = —2 and fj, = —^ then the point transformation 

x' = — 5-, t' = t, u' = — — , v' = - (3.15) 

v z 2x z u — v v 

maps the system (3.1) into itself. This latter point transformation forms a cyclic group of 
order 2. From this result we deduce that the integrated form of Eq.(2.6), 

\ 2 -i 

X \ IX 



v t = I — v xx + 

v x 2v x 



remains invariant under the point transformation x' = —x/v 2 , t' = t,v' = 1/v. 
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Case 3: If n / — 2, —1, and \x = ^w then the point transformation 

, g , , ?_ _J_ , J^(n + 2) 2n+2 

x'=t;2n + 2 f = t U = X n+2yn+l U L % l'=—± -L X n+2 (3.16) 

4(n+l) 2 
transforms the pdes 

/ / / ' \ ' /—(n+2) i x /-(ra+1) /„ t-n 

v x , = x u , v t i = X\x u u x i + X2U (3-17) 

2 

to the pdes (3.1), where Ai = if-xyya an d A2 = — ^f+rp • ^ we use ^ ne integrated forms 
of Eq.(2.6) and Eq.(2.2) we deduce that the point transformation 

, n n(n + 2) 2n+2 
X ' = V2n + 2 ; t =t, V = — -77X ™+ 2 

4(n + 1)^ 
maps the pde 

-(n+2) / / v. -(n+l) 



<> = Ai (f J ^ + (A2 - A:) (^ J 



into the pde 

n /v r \ n + 1 



4 More potential transformations 

We consider the standard nonlinear diffusion equation 

u' t , = (/(u'K')x' (4.1) 

which can be written, with the introduction of the potential variable v, as a system of two 
pdes 

v' x ,=u', v' t , = f(u')u' x ,. (4.2) 

Here we present point transformations of the form (3.9) which map Eqs.(4.2) into Eqs.(3.1) 
(or into Eqs.(3.2)). 

Case 1 : If n = — 1 the point transformation 

x' = v + 2t, t' = t, u ' = -^t + 1, v' = v + 2t + clnx, (4.3) 

ux z 

transforms the system (4.2) with / = —rz\, c a nonzero constant, into the system (3.2) 
where n = 1. 

Case 2: If n 7^ —2, —1 the point transformation 

, , , ?_ _-, n + 2 2(u+i) 

X = U, t = t, U =X n + 2 U , V = — r-X ™+ 2 , (4-4) 

2(n + 1) 

maps the system (4.2) with / = v! into the system (3.1) where u = ™+2 • ^ we 

use the integrated forms of Eq.(2.6) and Eq.(4.1) we deduce that the point transformation 

„ 2n + 2 

x' = v , t' = t, v' = 2n+2 x n+2 ma P s the pde 

/ _ / -(n+2) / 
v t' — v x' v x'x' 
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into the pde 



v x \ n n (v x \ n + 1 



v t = — I v xx . 

x J n + 2 \ x 

We note that if n = we obtain the hodograph transformation x' = v, t' = t, v' = x 
which connects the pde v' t , = v' x ,~ v' x , x , and the linear heat equation vt = v xx . This latter 
transformation is also a special case of the known general result where this later hodograph 
transformation connects the equations Vt = F{v x )v xx and Vt = (v x )~ 2 F(l/v x )v xx [9]. 
Furthermore point transformation (4.4) leads to the contact transformation 

dx' = uxdx + (xu) x dt, dt' = dt, v! = i~"+2u _ (4-5) 

which connects the nonlinear diffusion equation u' t , = I u'~ (n 'u x , ) and Eq.(2.6), where 
[i = +2 ■ W e observe that if n = 0, transformation (4.5) maps the well known nonlinear 
diffusion equation u' t , = (u'~ u' x , J into the linear pde ut = u xx + \u x . 
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